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Abstract 

When bifunctional transcription factors activate and repress target genes within 
the same cell, these opposing activities must be encoded in regulatory DNA. Here, we 
use cellular resolution gene expression data and computational modeling to investigate 
Hunchback (Hb) bifunctionality in Drosophila embryogenesis. Previous computational 
models predicted that Hb both activated and repressed the enhancer controlling even- 
skipped (eve) stripes 3 and 7 (eve3+7). We tested this hypothesis by measuring and 
modeling eve expression under multiple genetic perturbations and found that the 
eve3+7 enhancer could not explain endogenous stripe 7 behavior. To explain this 
discrepancy, we measured the response of an extended eve stripe 2 enhancer that 
drives expression of eve stripes 2 and 7 (eve2+7). We found that the behavior of 
endogenous stripe 7 is explained by the combined behavior of both enhancers, eve3+7 
and eve2+7. Bifunctionality arises from Hb activating the eve2+7 enhancer and 
repressing the eve3+7 enhancer. This pair can thus be considered "shadow enhancers" 
that both direct eve stripe 7, but respond to Hb in opposite ways. This example may 
illustrate a general way of encoding bifunctional regulation in the genome. 
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Introduction 

Transcription factors (TFs) are typically categorized as activators or repressors, 
but many TFs can act bifunctionally by both activating and repressing expression of 
their target genes (1-5). In cases where TFs activate and repress targets in the same 
cells, bifunctionality must be locally encoded in regulatory DNA sequence. Determining 
the regulatory DNA sequence features that control TF bifunctionality will advance two 
fundamental challenges of decoding transcriptional networks: predicting expression 
patterns from regulatory sequence, and deciphering how network topology dictates 
systems-level properties including gene expression precision and robustness to genetic 
and environmental perturbations. 

Here, we investigate how TF bifunctionality is encoded in regulatory DNA using a 
classic example: the Drosophila segmentation gene, hunchback (hb) (1 , 6-10). Hb both 
activates and represses the seven-striped even-skipped (eve) gene expression pattern 
by acting on multiple enhancers, genomic regions responsible for tissue specific gene 
expression (11). Hb activates eve stripes 1 and 2 and represses stripes 3,4,5,6 and 7 
(6, 8, 12). With our collaborators, we recently developed models of eve stripe regulation 
that suggested, consistent with previous models, that Hb bifunctionally regulates 
expression of eve stripes 3 and 7 (Fig. 1 ) (1 3, 1 4). Here, we test this hypothesis further 
using quantitative expression data in genetically perturbed embryos. We focus on 
whether Hb acts as both an activator and a repressor in the annotated eve3+7 
enhancer. We measured both the expression driven by the eve3+7 enhancer in a 
reporter construct and the endogenous eve expression pattern at high resolution in 
embryos with perturbed Hb expression levels. We then used these data to challenge 
our computational models. 

We found that Hb bifunctionality is encoded by separate enhancers that both 
direct eve stripe 7 expression. The first is the annotated eve3+7 enhancer, where Hb 
acts as a repressor. The second is an extended piece of regulatory DNA encompassing 
the minimal eve stripe 2 (eve2) enhancer that drives expression of eve stripes 2 and 7 
(eve2+7) (6, 15, 16). In the eve2+7 enhancer, Hb acts as an activator. Therefore, eve 
stripe 7 is controlled by a pair of shadow enhancers, separate sequences in a locus that 
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drive overlapping spatiotemporal patterns (17). Notably, these shadow enhancers 
respond to Hb in opposite ways and therefore use different regulatory logic. The 
separation of activation and repression into distinct enhancers may be a general 
mechanism of encoding TF bifunctionality in the genome. 
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Results 

eve enhancer reporter patterns do not match the endogenous eve pattern 

To determine if Hb bifunctionality is encoded in the annotated eve3+7 enhancer, 
we compared the pattern of the endogenous eve stripes to the pattern driven by a IacZ 
reporter construct in two genetic backgrounds. We refer to the resulting patterns 
throughout the manuscript as "the eve3+7 reporter pattern" and "the endogenous 
pattern" (Fig. 2). We examined both WT embryos and embryos where expression of Hb 
had been perturbed by removing bed, one of its key regulators, using RNAi (bed RNAi 
embryos) (18). We quantitatively measured expression patterns at cellular resolution 
using in situ hybridization, 2-photon microscopy, and an image processing toolkit 
developed specifically for Drosophila embryos (methods) (19, 20). We then averaged 
these data together into gene expression atlases (21). Importantly, the reporter 
construct isolates the activity of the annotated eve3+7 enhancer while the endogenous 
pattern integrates the activity of the whole locus. 

Our high resolution measurements revealed discrepancies between the 
endogenous pattern and the eve3+7 reporter pattern (Fig. 2). In WT embryos, the 
eve3+7 reporter pattern overlaps the endogenous eve stripes, but these stripes are 
broader, have uneven levels, and the peaks lie posterior to the endogenous peaks (Fig. 
2). These discrepancies were more pronounced in bed RNAi embryos than in WT 
embryos, especially for the anterior stripe (Fig. 2). When we tested reporters for other 
eve enhancers, we also found that they did not fully recapitulate the endogenous 
pattern (Fig. S1, S2). 

To test if the discrepancies between the eve3+7 reporter pattern and the 
endogenous pattern resulted from differences in eve and iacZ mRNA half-lives, we 
measured the expression driven by an eve locus BAC reporter where the coding 
sequence had been replaced with IacZ (a generous gift from Miki Fujioka). In both WT 
and bed RNAi embryos, the peak positions and widths of the BAC reporter pattern were 
more faithful to the endogenous eve pattern, but still did not match exactly (Fig. 2, S1, 
S2). Differences between the endogenous and BAC reporter patterns must arise from 
differences in the transcripts. Differences between the BAC reporter and the isolated 
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enhancer reporter patterns must arise from regulatory DNA not included in the reporter 
constructs. These data indicate that the eve3+7 reporter construct is missing relevant 
regulatory DNA. Additional regulatory DNA in the endogenous locus may respond to 
other TFs or respond to the same TFs differently leading to differences in how the 
endogenous pattern is computed from the concentrations of its regulators. We therefore 
hypothesized that the endogenous locus and the eve3+7 enhancer perform different 
computations to produce eve stripes 3 and 7. We next tested this hypothesis using 
computational models of eve regulation. 

Computational models suggest that Hb activates and represses endogenous eve 
stripes 3 and 7, but only represses the eve3+7 enhancer 

With our collaborators, we previously identified two empirical computational 
models for the expression of eve stripes 3 and 7. (Fig. 1) (13). These models used 
logistic regression to directly relate the concentrations of input regulators to output 
expression in single cells by fitting a parameter for each regulator that reflects both TF 
strength and TF-DNA interactions (13). In the linear model, Hb has one parameter and 
only represses. In the quadratic model, Hb has two parameters, and activates at low 
concentrations while repressing at high concentrations (13). Both models perform 
equally well in WT embryos, but they make different predictions under genetic 
perturbation. Specifically, the quadratic model predicts published data from hb ventral 
misexpression experiments (see below). 

To investigate patterns driven by individual enhancers, these models were fit on 
data parsed from the entire endogenous eve pattern because data for enhancer 
reporters was not available (13). To relate the results to individual enhancers, we 
employed a standard assumption: the endogenous expression of eve stripes 3 and 7 
could be attributed to the activity of the annotated eve3+7 enhancer. Here, we test this 
assumption explicitly by measuring and modeling the eve3+7 reporter and endogenous 
patterns separately. 

We compared the performance of the linear and quadratic models in WT and bed 
RNAi embryos. As input regulators we used Hb protein and gt, til and kni mRNA, and 
we used thresholded endogenous or eve3+7 reporter mRNA data as our target output 
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patterns (Fig 1, methods). All of these regulators, especially Hb, are perturbed in bed 
RNAi embryos (Fig. S3) (18). We report our modeling of the third time point, which is 
representative of results for other time points (Fig. S6), and evaluated model 
performance by computing the area under the receiver operating characteristic curve 
(AUC) (22). 

We first addressed the endogenous eve pattern: we fit our models in WT 
embryos and used the resulting parameters to predict expression in bed RNAi embryos. 
In this case, we found that the quadratic model more precisely predicted the perturbed 
endogenous eve pattern. Both models correctly predicted the positional shifts of stripe 7 
and a wide anterior stripe, but the quadratic model was more accurate than the linear 
model (AUCiinear = 0.93, AUCquad = 0.98) (Fig. 3F, Fig. S4). These analyses indicated 
that the quadratic model captured the activity of the whole locus by allowing Hb to both 
activate and repress eve stripes 3 and 7. 

We next addressed the eve3+7 reporter pattern: again, we fit our models in WT 
embryos and used the resulting parameters to predict expression in bed RNAi embryos. 
To our surprise, the linear model was more accurate than the quadratic model in this 
case (AUCiinear = 0.90, AUCquad = 0.87) (Fig. 3L). Although neither model captured the 
dorsal-ventral modulation of the pattern, the linear model accurately predicted the 
posterior boundary of the anterior stripe. We controlled for several factors that may 
confound model performance. We assessed sensitivity to changes in regulator 
concentrations, refit the models in bed RNAi embryos alone, and refit the models on all 
of the data, none of which changed our conclusions (Fig. S5 and S6, Supplemental 
Note 1). Although these differences in model performance were subtle, the results 
supported our hypothesis that the eve3+7 reporter pattern is regulated differently than 
the endogenous pattern. Specifically, they suggested that the endogenous pattern 
required bifunctional Hb while the eve3+7 reporter pattern required only Hb repression. 

Based on these results, we hypothesized that Hb activation is encoded in 
regulatory DNA outside the annotated eve3+7 enhancer. The differences in model 
performance were not conclusive of their own accord, prompting us to seek another 
perturbation to validate this hypothesis. We therefore returned to the perturbation that 
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previously distinguished the linear and quadratic models, ventral mis-expression of hb 
(13, 14). 

hb mis-expression confirms that the endogenous eve pattern and the eve3+7 
reporter pattern respond to Hb differently 

In llsley et al., we preferred the quadratic model because it qualitatively predicted 
the behavior of a classic mis-expression experiment. Mis-expressing hb along the 
ventral surface of the embryo (sna::hb embryos) causes eve stripe 3 to retreat and bend 
and stripe 7 to bend and expand (23). In simulations of this perturbation, the quadratic 
model predicted this behavior while the linear model did not (Fig. 4 E and F reproduced 
from (13)). We hypothesized that if Hb bifunctionality is encoded in the whole locus but 
not in the annotated eve3+7 enhancer, the endogenous and eve3+7 reporter patterns 
would respond differently to hb misexpression. To test this hypothesis, we repeated this 
perturbation and measured both the endogenous eve pattern and the eve3+7 reporter 
pattern quantitatively at cellular resolution (Fig. 4). 

In sna::hb embryos, the endogenous pattern behaved differently from the 
eve3+7 reporter pattern. As previously observed, the endogenous eve stripe 3 retreated 
from the ventral Hb domain and bent posteriorly; the endogenous stripe 7 did not retreat 
at all, and instead expanded and bent anteriorly. Both of these behaviors were predicted 
by the quadratic model applied to simulated mis-expression data (Fig. 4 B and E). By 
contrast, in the eve3+7 reporter pattern, both stripes 3 and 7 retreated from the ventral 
Hb domain, consistent with the predictions of the linear model applied to simulated mis- 
expression data (Fig. 4 C and F) (13). Thus, under this mis-expression perturbation, the 
linear model predicted the behavior of the eve3+7 reporter pattern while the quadratic 
model predicted the behavior of the endogenous pattern. The subtle quantitative 
difference between the two models that we saw in bed RNAi embryos was corroborated 
by a strong qualitative difference in sna::hb embryos. These data indicate that Hb 
represses the eve3+7 reporter pattern and both activates and represses the 
endogenous pattern. Hb repression is encoded in the eve3+7 enhancer; we next 
sought to identify the regulatory DNA in the eve locus that directs stripe 7 expression, 
but is activated by Hb. 
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Two shadow enhancers encode bifunctional Hb regulation of eve stripe 7 

We hypothesized that Hb activation of eve stripe 7 was encoded in regulatory 
DNA near the eve2 enhancer. We focused on the eve2 enhancer for several reasons: 
Hb is known to activate the eve2 enhancer (7); longer versions of eve2 drive stripe 7 in 
some embryos (15, 16, 24, 25); orthologous eve2 enhancers from other species 
sometimes drive stripe 7 expression (26, 27); finally, in sna::hb embryos the border of 
the expanded stripe 7 appeared to be set by Kruppel (Kr), a known regulator of eve2 
(Fig. S7) (6, 28). We measured expression driven by an extended version of the 
minimal eve2 enhancer construct that drove a robust stripe 7 pattern (Table S2); we call 
this enhancer reporter construct eve2+7. Since both the eve3+7 and eve2+7 enhancers 
drive stripe 7 expression, they can be considered a pair of shadow enhancers for stripe 
7(17). 

We found that Hb activates stripe 7 expression in the eve2+7 enhancer. In 
snar.hb embryos, the stripe 7 region of the eve2+7 reporter pattern expanded, 
recapitulating the bulging behavior observed in the endogenous eve pattern (Fig. 4B 
and D). We conclude that Hb activates endogenous eve stripe 7 through the eve2+7 
enhancer. Taken together, our results indicate that Hb controls eve stripe 7 expression 
by activating and repressing distinct enhancers. 
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Discussion 

We uncovered a pair of shadow enhancers in the eve locus that both direct 
expression of stripe 7 using distinct regulatory logic: one uses Hb as an activator, the 
other uses Hb as a repressor. We measured expression of the endogenous eve locus 
and transgenic reporter constructs at cellular resolution under two genetic perturbations. 
We contextualized our results by comparing two computational models with different 
roles for Hb: a linear model where Hb is a dedicated repressor and a quadratic model 
where Hb is a bifunctional regulator (13). Guided by the modeling, we found that stripe 7 
is encoded by two enhancers: Hb represses the eve3+7 enhancer and activates the 
eve2+7 enhancer. These two shadow enhancers therefore use Hb in different ways to 
position overlapping patterns. This form of regulatory redundancy may be a general way 
to encode TF bifunctionality. 

Expression patterns driven by reporter constructs do not precisely match the 
endogenous pattern 

"Veteran enhancer-bashers, and those who carefully read the papers, know that 
'minimal' enhancer fragments do not always perfectly replicate the precise spatial 
boundaries of expression of the native gene..." (17). Our data clearly support this often 
neglected aspect of enhancer reporter constructs. One explanation offered for such 
discrepancies is different mRNA half-lives. We controlled for this possibility with a BAC 
reporter where the eve coding sequence was replaced with lacZ and found better, but 
not perfect, agreement with the endogenous pattern. We conclude that transcript 
sequence features contribute to the differences between reporter and endogenous 
patterns, but that additional regulatory DNA in the locus also plays a role. This result 
highlights the limitations of enhancer reporter constructs for recapitulating behavior of 
endogenous loci and the importance of using BAC reporters or genomic editing to study 
loci with multiple enhancers (29-31). 
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eve stripe 7 is encoded by two shadow enhancers 

Early efforts to dissect the regulatory architecture of the eve locus suggested 
that stripe 7 activity was distributed over DNA encompassing both the eve3+7 and eve2 
enhancers (15, 16). We find that there are at least two regions of regulatory DNA that 
position stripe 7 using different regulatory logic. The classically annotated eve3+7 
enhancer is repressed by Hb (8, 23, 32), while the eve2+7 enhancer, which 
encompasses the minimal eve2 enhancer, is activated by Hb. We think that activation 
of eve2+7 is direct because it is clear from in vivo DNA binding data and binding site 
mutagenesis that Hb directly activates the minimal eve2 enhancer (7, 33). The 
redundancy in eve stripe 7 regulation may confer robustness to genetic or 
environmental stresses (34-36), may increase synchrony or precision (37), may 
facilitate temporal refinement of patterns (38), or may arise from genetic drift (39). 

It is likely that these two enhancers are differentially sensitive to additional TFs. 
Previous studies have revealed the eve3+7 enhancer is activated by the spatially 
uniform TFs Stat92E and Zelda, the anterior boundary of stripe 7 is set by Kni 
repression, and the posterior boundary is set by Hb repression (8, 23, 32, 40). The 
minimal eve2 enhancer is activated by Bed and Hb, its anterior boundary is set by Slp1 
and Gt, and its posterior boundary set by Kr (7, 24, 28, 41). Kr appears to set the 
boundary of both the expanded endogenous eve stripe 7 and the eve2+7 reporter 
pattern in sna::hb embryos (Fig. S7), but eve3+7 is not sensitive to Kr, as eve stripe 3 
sits directly beneath the Kr pattern and this enhancer has no predicted Kr binding sites 
(Fig. S7). In agreement with others, we speculate that the anterior boundary of eve 
stripe 7 in eve2+7 may be set by Gt (25). However, we cannot rule out the possibility 
that this boundary is set by limiting levels of Hb activation or by Kni repression. 

The molecular mechanism by which Hb represses eve3+7 and activates eve2+7 
remains unclear. One hypothesis is that other TFs convert Hb from a repressor into an 
activator. For example, there is experimental evidence for activator synergy between 
bed and hb (42) and Hb/Cad activator synergy has been proposed based on 
computational work (43). Another proposed mechanism is that monomeric Hb is an 
activator, but Hb dimers are repressors (14, 44). High and low Hb concentrations may 
also be correlated with some other spatially varying factor in the embryo, such as 
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phosphorylation by the MAPK pathway in the poles (45), but data from sna::hb embryos 
reduces the likelihood of this possibility. Testing these hypotheses will require 
quantitative data in additional genetic backgrounds and mutagenesis of individual 
binding sites in the two enhancers. 

The quadratic model is a superposition of two computations 

Models are not ends in and of themselves, but merely means to formalize 
assumptions and develop falsifiable hypotheses (46, 47). The quadratic model, which 
includes Hb as a bifunctional TF, accurately predicts the behavior of the locus in all cells 
of WT and perturbed embryos, but it does not predict the behavior of either individual 
enhancer. The interpretation in llsley et al. that Hb bifunctionality is a feature of the 
canonical eve3+7 enhancer was based on a common assumption: that the endogenous 
expression pattern could be attributed to the annotated enhancer. Here we show that 
Hb bifunctionality is encoded in separate enhancers. The quadratic model works 
because it combines the critical features of the eve3+7 enhancer and the eve2+7 
enhancer, effectively behaving as a superposition of the two activities. In the future, we 
plan to develop computational models of each enhancer and uncover how these two (or 
more) activities are combined. 

Conclusion 

We tested the hypothesis that Hb bifunctionality is encoded in the eve3+7 
enhancer and discovered that it is actually encoded in two separate enhancers that 
respond to Hb in opposite ways. We show that expression patterns driven by annotated 
enhancers differ from the endogenous pattern, especially under perturbation, and that 
these differences can be due to relevant, yet unannotated, regulatory DNA. The stripe 7 
shadow enhancers reside in a classic example of a modular locus (48), implying that 
regulatory complexity may be pervasive. Since the enhancers are active in the same 
cells, Hb bifunctionality must be encoded in their DNA sequences. This example 
provides an opportunity to uncover sequence features governing Hb bifunctionality, 
which will improve our ability to interpret regulatory DNA and infer connections in gene 
regulatory networks. 
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Materials and Methods 
Fly Work 

The bed RNAi gene expression atlas is described in Staller et al. 2014 
(submitted) and available at depace.med.harvard.edu. Briefly, we combined short 
hairpin RNA knockdown of bed with in situ hybridization and 2-photon imaging and 
automated image segmentation (21, 49-51). Hb protein stains used a guinea pig an\\-hb 
from John Reinitz (University of Chicago, IL). Embryos were partitioned into six time 
points using the degree of membrane invagination (0-3%, 4-8%, 9-25%, 26-50%, 
51-75%, and 76-100%) which evenly divide the ~60 min blastoderm stage (20). All 
enhancer reporters are in pBOY and integrated at attP2 (26, 52) (Table S1). The eve 
locus bacterial artificial chromosome (BAC) lacZ reporter was a gift from Miki Fujioka 
(Thomas Jefferson University, PA). It begins 6.4 kb upstream of the eve transcriptional 
start site (TSS) and ends 11.3 kb downstream of the eve TSS. The eve coding 
sequence has been replaced by lacZ and the adjacent gene, TER94, has been fused to 
GFR It is effectively a reporter for the whole eve locus, hb ventral misexpression was 
performed as described in Clyde et al., 2003 using two copies of the transgene on 
chromosome 2. 

Building the coarsely aligned snar.hb gene expression atlas. 

We determined the genotype of the snar.hb embryos by examining the eve or 
fushi-tarazu (ftz) mRNA patterns. Embryos were aligned morphologically to create a 
coarsely registered gene expression atlas (21). Data is available at 
depace.med.harvard.edu. 

Logistic modeling of enhancer gene regulatory functions 

The logistic modeling framework was developed and described in detail 
previously (13). All modeling was performed in MatLab (MathWorks, Natick, MA) using 
the DIP image toolbox (diplib.org) and the PointCloudToolBox (bdtnp.lbl.gov). Ilsley et 
al. used protein data for Hb and Gt, whereas we used Hb protein and gt mRNA data. 
For genes where we used mRNA data, the mRNA and protein patterns are known to be 
correlated (21 , 53). For the enhancer lacZ reporters, we thresholded cells to be ON or 
OFF by creating a histogram of the expression data (50 bins), identifying the bin with 
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the most counts and adding one standard deviation. Our ON set included all cells 
expressing the reporter, and our OFF set includes all other cells. All regulators are 
maintained as continuous values. 

To threshold the endogenous WT eve pattern into ON and OFF cells we used 0.2 
for all time points (13). To threshold the endogenous eve patterns in the bed RNAi atlas, 
we used the lowest threshold that would separate the stripes: 0.1, 0.15, 0.15, 0.2, and 
0.21 for T=2 through T=6 respectively. To compare the modeling of the reporter and the 
endogenous patterns, the ON set included all cells in the endogenous eve stripes 3 and 
7 and the OFF set included all other cells. This OFF set is different from llsley et al., but 
this change does not have a large effect on the model prediction AUC scores in bed 
RNAi embryos (Table S1). 
Sensitivity analysis 

For the sensitivity analysis (Fig. S5), for each TF, we scaled the concentration of 
the bed RNAi atlas in silico and recomputed the model AUC scores. 
Binding site predictions 

For the Kr binding site analysis in Fig. S7, we predicted binding sites using 
PATSER (stormo.wustl.edu) with a position weight matrix derived from bacterial 1 -hybrid 
data (54). Binding sites were visualized using InSite (cs.utah.edu/~miriah/projects). 
Quantifying concordance between reporters and endogenous patterns 

For each embryo, we used the pointcloud toolbox in Matlab to find pattern 
boundaries by creating 16 anterior-posterior line traces and finding the inflection point of 
each. Finding the boundary by using half the maximum value of the stripe peak 
identifies a very similar boundary to the inflection point. To find the peaks of the 
endogenous and reporter stripes, we took one line trace along the lateral part of the 
embryo using the pointcloud toolbox and found the local maxima. 
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Figure Legends 

Figure 1 : 

The linear and quadratic models formalize two alternative regulator sets for eve stripes 
3 and 7. (A) The linear model includes repression (red) by Hb, knirps (kni), giant (gt), 
and tailless (til) and activation (blue) by a constant term that represents spatially uniform 
factors. The quadratic model includes activation by a linear Hb term and repression by 
a quadratic Hb term, kni, til, and uniform factors. (B) A schematic of the logistic 
regression framework. Logistic regression calculates the probability the target will be 
ON based on a linear combination of the concentrations of regulators (/v). We fit models 
in WT and use the perturbed regulator gene expression patterns to predict the 
perturbed eve patterns in bed RNAi embryos. 

Figure 2: 

The eve3+7 reporter pattern differs from the endogenous pattern. (A) The eve locus 
contains 5 annotated primary stripe enhancers. The endogenous pattern integrates the 
activity of the whole locus. The BAC reporter construct also integrates the activity of the 
whole locus, but the transcript is the same as the eve3+7 reporter construct. The 
eve3+7 reporter construct isolates the activity of the annotated enhancer sequence. (B) 
WT expression patterns are represented as line traces where anterior-posterior (A-P) 
position is plotted on the X-axis with expression level on the Y-axis for a lateral strip of 
the embryo. Endogenous eve pattern (gray), eve3+7 reporter pattern (red). The reporter 
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pattern was manually scaled to match the level of the endogenous pattern. (C) Line 
traces in bed RNAi embryos. (D) The boundaries of the endogenous pattern (gray), the 
eve3+7 reporter pattern (red), and the BAC reporter pattern (blue) at T=3. All error bars 
are the standard error of the mean. The BAC reporter pattern is more faithful to the 
endogenous pattern than the eve3+7 reporter pattern, especially in the anterior of bed 
RNAi embryos (eve 3/7 ant). The endogenous pattern is shaded for visual clarity. (E) 
Peak positions of stripes 3 and 7, calculated from the line traces in B and C. The 
eve3+7 reporter pattern shows better agreement to the endogenous pattern in WT than 
in bed RNAi embryos. (F) Stripe widths, calculated from the inflection point of the line 
traces in B and C. The eve3+7 reporter pattern is wider than the corresponding 
endogenous pattern. 

Figure 3: 

In bed RNAi embryos, the quadratic model more accurately predicts the endogenous 
pattern, and the linear model more accurately predicts the eve3+7 reporter pattern. (A) 
The endogenous eve pattern in WT embryos. Cells with expression below an ON/OFF 
threshold (methods) are plotted in gray. For cells above this threshold, darker color 
indicates higher level. (B) The predictions of the linear and quadratic models in WT 
embryos. (C) Comparison of model predictions to the endogenous pattern in WT 
embryos. Green cells are true positives, purple cells are false positives, dark gray cells 
are false negatives, and light gray cells are true negatives. For visualization, the 
threshold is set to 80% sensitivity, but the AUC metric quantifies performance over all 
thresholds. (D) The endogenous eve pattern in bed RNAi embryos. (E) The predictions 
of the linear (L) and quadratic (Q) models in bed RNAi embryos. (F) Comparison of 
model predictions to the endogenous pattern in bed RNAi embryos. The quadratic 
model more accurately predicts the endogenous pattern in bed RNAi embryos. (G-L) 
Same as A-F, respectively, for the eve3+7 reporter pattern. The linear model predicts 
the eve3+7 reporter pattern more accurately in bed RNAi embryos. Model parameters 
are in Table S1 . 
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Figure 4: 

In hb ventral misexpression (sna::hb) embryos, the quadratic model predicts the 
endogenous pattern while the linear model predicts the eve3+7 reporter pattern. (A) The 
distribution of Hb in WT and sna::hb embryos, from a lateral view (left) and ventral view 
(right). The expression level of each regulator is shown for individual cells: cells with 
expression below an ON/OFF threshold (methods) are plotted in gray. For cells above 
this threshold, darker colors indicate higher levels. (B) Endogenous eve pattern. (C) 
The eve3+7 reporter pattern. Both stripes retreat from ectopic Hb. (D) The eve2+7 
reporter pattern. Stripe 7 expands into the ectopic Hb domain. (E-F) Bottom (ventral) 
view of predictions of the quadratic model (E) and linear (F) models based on simulated 
snar.hb data. OFF cells are light pink and ON cells are red. Reproduced from llsley et 
al. 2013. 
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